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Abstract. The Andreev reflection probability for a ferromag- 
net /isolator /superconductor (FIS) contact at the arbitrary spin- 
dependent amplitudes of the electron waves transmitted through 
and reflected from the potential barrier is found. It is shown that 
Andreev reflection probabilities of electron and hole excitations in 
the FIS contact are different. The energy levels of Andreev bound 
states are found. The ballistic conductance of the point FIS contact 
is calculated. 

PACS: 74.50.+r, 74.80.-g, 75.30.Et 

One of the manifestations of the exchange field in a ferromagnetic 
metal (F) is the presence of electron spin subbands with different 
values of Fermi momenta: p-j- for the subband with spin up and 

for the subband with spin down. As a consequence, for layered 
F/S structures (S stands for superconductor) the spatial depen- 
dence of the anomalous Green's function (GF) in a ferromagnet has 
an oscillatory character. One of the impressive manifestations of 
such oscillations and related phase shifts is a recent observation of 
spontaneous zero-field supercurrents at temperature lower than the 
junction 0—tt transition temperature in superconducting networks of 
SFS junctions with weakly ferromagnetic barriers pQ. The influence 
of the oscillatory character of the anomalous GF in a ferromagnet 
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on the properties of various hybrid F/S structures is studied well 
enough (see reviews [2-4]). 

Another consequence is the suppression of Andreev reflection [5] . 
When a polarized electron from the subband with, for example, spin 
up gets into a superconductor, the reflected hole moves into the 
subband with spin down. Consequently, the efficiency of Andreev 
reflection is determined by the number of conducting channels in 
a subband with a smaller value of the Fermi momentum. As a 
result, the subgap conductance of an F/S contact decreases with 
the increase of the polarization of a ferromagnet j6]. 

Effects of spin filtering [TJ, [8], [9] and spin mixing [10] are man- 
ifested in the dependence of moduli and phase shifts of the am- 
plitudes of electron states on the Fermi surface reflected from r a 
(r a = ^/R a exp(i9^) ; R a = 1 — D a ) and transmitted through a 
potential barrier d a (d a = \/D a exp (i #„)) on a (a = f,| is the 
spin index). These effects are the consequence of the presence of 
the exchange field in a ferromagnet as well. 

The possibility to study the influence of the spin mixing effect 
on the I — V characteristics of superconducting weak links con- 
taining a magnetically active interface appeared after the boundary 
conditions (BCs) for the quasiclassical GF were obtained. 

In paper BCs for the quasiclassical GF for two metals in 
contact via a magnetically active interface in terms of an interface 
scattering matrix were derived. These equations were solved for a 
junction in the tunneling limit [11] and for a contact of a supercon- 
ductor with a ferromagnetic insulator |12j . In paper [TO] , BCs for the 
retarded and advanced quasiclassical GFs were obtained in terms of 
Riccati amplitudes |13j . [14] . In paper [15], BCs in terms of Riccati 
amplitudes were obtained for the nonequilibrium quasiclassical GF. 

The equations, obtained in papers [JO] and [T5], were solved for 
magnetically active interfaces with finite transmission (for SFS [5], 
HO], for NFS [T5] (N stands for normal metal), for S-FIF-S [IS]). 
These solutions show that Andreev bound states appear within the 
superconducting gap [8] , [10] , [15] , and the — it transition in the 
SFS junction is possible [8], [TO] . 

In papers [9] and [TTJ, quasiclassical equations of superconduc- 
tivity for metals with a spin-split conduction band were derived and 
BCs for the temperature quasiclassical GF for the F/S interface were 
obtained. The model interface was the same as in [TT], [18] . 
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The aim of this work is to study the influence of spin-dependent 
phases of the amplitudes of the electron states reflected from and 
transmitted through a potential barrier on Andreev reflection in a 
point FIS contact. 

Calculations are carried out by the method of quasiclassical GFs 
with BCs for GFs obtained in papers [3], [17]. Below the dependence 
of the Andreev reflection probability on spin-dependent phase shifts 
9 ^ and 9 r a will be found and the results of the numerical calculation 
of the dependence Gfis(V) for a rectangular potential barrier and 
ferromagnets with high polarization will be discussed. 

1. Differential conductance of a point FIS contact. In 
various hybrid F/S structures Andreev reflection is modified. The 
reflected hole has some parameters (for example, the velocity mod- 
ulus and phase shift) different from those of the incident electron 
because it moves in a subband with the opposite spin. Such spin- 
discriminating processes due to the exchange field in a ferromag- 
net lead to the formation of Andreev bound states inside the gap 

mm- 

The enegy of Andreev bound states depends on the spin index 
[8], [10]. As a result, the spectral density of condauctance Gfis of 
the FIS contact at zero voltage is no longer a symmetrical function 
of energy e. The condition of the time reversal invariance has the 
form Gpis( £ , a ) — Gpis(~ s i — ex)- The generalization of the con- 
ductance Gpisiy) 0) [IH] for this case results in the following 
expression for Gpisiy)'- 
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In Eq. (CD), A is the contact area; f z is the Pauli matrix; pn 
is the momentum in the contact plane; (g s , T s ) and (g a , T a ) are 
quasiclassical retarded (R) and advanced (A) GFs symmetric and 
antisymmetric with respect to the projection of the momentum p 
on the Fermi surface on the axis x, respectively [9]. Calculations 
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in Eq. ([I]) are to be carried out on the boundary of any contacting 
metal. 

2. Finding GFs and conductance. Let us assume that the 
barrier with the width d is located in the region a < x < b {d — 
b — a), the superconductor occupies the region x > b, and the 
ferromagnet occupies the region x < a. To find GFs, for each metal 
one has to solve quasiclassical equations of superconductivity for 
metals with a spin-split conductivity band simultaneously with their 
BCs derived in paper [9]: 

did 
sign^.)— # + - vy — (v~ g + gv~ ) + [K, g}_ = 0, 
ox 2 op 

K = -iv x 2 (ie n T z + A - T)v x 2 - i(p x - r x p x i x )/2, 

[a, 6]_ = ab — 6a. (2) 

In this section, e n = (2n + l)nT is the Matsubara frequency; X 
is the self-energy part; g are matrix temperature GFs: 
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where A is the order parameter, and p x is the projection of the 
momentum on the Fermi surface on the axis x. Matrices v have the 
same structure as p x . 

BCs for the specular reflection of electrons from the boundary: 
p\\ = p[ sini?| = p] sini?f = ps sin ■&$■, have the form [9]: 

(l!) d = (If )d, (^ah = (^aU 



Ra-V R- a )(T+) n = a 3 (g a ) n , 



7 = v R a (gt)d + «i(^)n, 
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- f + = (R a ) 5 Q s ) d + a 2 Q s ) n , (3) 

where g± {s) = 1/2 [<? o(s) ± g a{s) ]. Functions T o(s) are determined 
analogously. The index d denotes the diagonal and n the nondi- 
agonal part of the matrix T d ^ = 1/2 [T ± t z Tt z ]. Coefficients a, 
are: 

i + v^rIt/d^dI 

«1(2) = 



a 3 (4) = 1 - v 7 ^^! ± V^Dl )• 

One can exclude GFs T ^ and T f from these relations and obtain 
a system of BCs only for the GF g [T7] : 



#a & 1 + Ma + #a & 3 + Ma = &3 ~ &4, 

g a h + 6 2 ^ a + l a h + &4<7 a =bi-b 2 . 



(4) 



Matrices fej in Eq. (jlj) are: 



h = T+# s + f s g s , 6 2 = # s f s + # s f - 
h = T+|+ + T s g~, b A = g + s T+ + g s f~ 6 



(5) 



GFs g are connected with GFs being solutions of Eq. ([2]) by the 
following relationships [5]: 
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0a=^ Y ^-(9 d a -9l a ); 0;=-i_z2. (6) 

The explicit form of functions T is not needed. These functions are 
found from BCs. The diagonal parts of matrices g are equal to the 
corresponding matrices g. Equations (T5]) for the ballistic contact are 
solved in paper [9], [18]. At the boundaries for x = b and x = a we 
have: 

g s s = go + go £?; gf = go - go €■ (?) 
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Matrices go are values of GFs g away from the boundary: 



9rJ 



9 S o 



9oT z + 



sign(£ n J r z 
1 



(8) 



-iA 



After the substitution of functions gf and <?f , expressed via g% and 
<?f by Eq. (J7j), in the system of BCs Eq. (HI) and their solution in the 
linear approximation with respect to the functions g% g£, we find 
the function g£: 
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Z=(l- y/R^) [9o ^s(6 a ) + i sin(0 Q )] (9) 
+ (1 + y/R^Ri) sign(e n ) [cos(9 a ) +ig$ sin(0 a )]. 

From Eqs. ([3]) and (jlj) we find the rest functions necessary to cal- 
culate conductance Eq. ([T]) and calculate conductance at the ferro- 
magnet side. 

After carrying out the analytical continuation in these functions 
(substitution % e n for e ± i5 for retarded and advanced GFs, respec- 
tively), we obtain the expression for the conductance cr F /s(V): 
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W = y/R^Ri; £ = v^ 2 - |A| 2 - (10) 

At 9 a = the expression for conductance obtained in paper [5] 
follows from Eq. (1101) . In the case of nonmagnetic metal, when 
= this expression is the same as that obtained in paper [IB] , 
and for D = 1/(1 + Z 2 ) this expression is the same as that obtained 
in paper [T9] . 

3. Andreev reflection. The quasiclassical GFs entering Eq. 
(jHJ) enable the conclusion that 

[1 " 9 s r z gff z -gtf z g^f z + t ff z f f t z 

-tjf,fff,] = 4[-^fj?f,]~i. (ii) 

Now, the comparison of the form of under-gap conductances in 
Eq. ([I]) and that of the corresponding Eq. (25) in paper [T3] shows 
that the matrix elements (g^) F an d iJ)a) F are the amplitudes of the 
Andreev reflection probability a(e, 9 a ) in FIS contacts for energies 
less than |A| (e 2 < |A| 2 ). Let us assume that o(e, 9 a ) are matrix 
elements of {g^Y- 



a(e, 9 a ) = V J z l = a(e, 9 a )e~ 1 ^ (12) 

Z= (1 - yjR<\R$\e cos(# Q ) - a/|A| 2 - £ 2 sin(0 Q )] 
+i (1 + v/ip^)[v/|A| 2 -£ 2 cos(# Q ) + e sin(0 a )]. 

The presence of the imaginary part in functions a(e, 9 a ) means that 
Andreev reflection is accompanied by the phase shift. The Andreev 
reflection probability A(e, 9 a ) (A(e, 9 a ) = a(s, 9 a )a*(e, 9 a )) is: 

A(e, 9 a ) = (13) 

Z = [1- v^p!] 2 |A| 2 

+ 4 v 7 ^! M A I 2 - cos( 9 a ) + e sin{9 a )} 2 . 

It follows from this equation that: (1) in terms of paper [10J spin- 
mixing angle for FIS contact is equal to 9 a (for SFS and NFS 
contacts Q = 9\ - 9\ = 0f - 0f |SJ , [EH , [E] ) ; (2) for 9 a < the 
Andreev reflection probability of the electron excitation with the 
spin projection a is larger than that of the hole excitation; for 9 a > 
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the Andreev reflection probability of the hole excitation with the 
spin projection a is larger than that of the electron excitation; (3) 
the Andreev reflection probability has maxima at e = (at the 
values of the energy of electron (hole) excitations corresponding to 
the energy levels of Andreev surface bound states) 

f e = |A| cos(# Q ) for 9 a < 0, 
€b ~ \e = -|A| cos(# Q ) for 6 a > 0. 1 j 

Below the results of the numerical calculations of phase shifts and 




Figure 1: Dependence of the phase shifts of reflection and transmission ampli- 
tudes on cos(??i). Lines with numbers 1,2, and 3 depict these dependences on 
cos(tf x ): T ; (0f -Of) and (6^ - 9 J), respectively. 

conductance are presented. In the numerical calculations the rela- 
tion between Fermi momenta of contacting metals was the following: 
Ps = (pt + Calculations are carried out for a rectangular 

barrier with the height C/ counted off the bottom of the conduc- 
tion band of a superconductor; [x(x) is the wave function of an 
electron in an isolator, xi x ) = Ci exp(7x) + C2 exp(— 7a;); 7 = 
^//c 2 + p|; A; 2 = 2m b (U — Ep) ; is the Fermi energy of a super- 
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Figure 2 : Dependence of the normalized conductance up/s(V)/ao from Eq. (|10p 
on the applied voltage for different values of the polarization of a ferromagnet 
S = p l /p 1 at the ratio A d (T)/2T =6. 



conductor, m;, is the mass of an electron in a barrier]. In this case 
the expressions for 9 d a and 9 r a have the following form: 

^ = eJ + i(pJ a o-pf6); 9 r a = 9 r a + 2ipF a a 

nd + f (Px, a Px-7 2 )tanh( 7 ci) ^ 

6> Q = arctan — ■ ^ (15) 



9L — arctan 



'^lV F x Jl 2 + (P S X ) 2 } tanh(7c/) , 



Z 

Z = l 2 M) 2 - «J 2 ] + [7 2 -pf < Q ] 2 tanh 2 ( 7 rf), 

so that the angle 9 a [9 a = {0 r a - 91 J/2 - - 9 d a )] = (9 r a - 
9^ a )/2 — (9 d — 9 d a ) does not depend on the location of the barrier. 

Figure 1 shows the dependences of the phase shifts on cos($j,). 
All angles are connected by specular reflection pii = pjsini!^ = 
PT sin i?f = ps sin $5. 
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The phase shift 9<[ slowly decreases as the polarization of the 
ferromagnet 5 decreases [from (-1.3) at 5 = 0.05 to (-1.5) at 5=0.5 
(p T d = 1; k/ Pl = 0.8) ] and [from (-0.7) at 5=0.05 to (-1.2) at 
5 = 0.5 (pid = 1; k/pi = 0.2 )]. It means that the points e& ap- 
proach zero as the polarization decreases and k/p^ increases, how- 
ever, at k/pi > 1, 0f rapidly decreases down to zero. 6. With 
the increasing parameter p^d and other parameters fixed (but for 
k/pi < 1) the angle 9 a also tends ton/ 2. Note that the spin-mixing 
angle 9 a for ferromagnets with large polarization is practically the 
same for all electron trajectories. 

The upper panel in Fig. 2 shows the results of the numerical cal- 
culations carried out according to Eq. ffTOl) not taking into account 
(dashed lines) and taking into account (solid lines) the phase shift 
9 a. The peaks in the dependence of the conductance on V (Fig. 
2, upper panel) correspond to the motion of the energy levels of 
Andreev surface bound states towards each other as the parameter 
fc/p-f increases. 

The lower panel in Fig. 2 shows the suppression of Andreev 
reflection due to the reduction of the number of conducting channels 
in a subband with a lower value of the Fermi momentum and the 
effect of spin filtering. Andreev surface bound states are formed in 



Figure 3: Structure of the diagrams corresponding to Andreev reflection in the 
superconductor: diagram a) one-act process; diagram b) two-act process. The 
vertex O is Andreev reflection of electronlikc (solid lines) and holelikc (broken 
lines) quasipaticles by the pair patential. The vertex • is the normal reflection 
of electronlike and holelike quasipaticles by the barrier potential. When the 
solid line transforms into the broken line, O denotes the vertex /3„ _ a . When 
the broken line transforms into the solid line, Q denotes the vertex /3_£ a . 
Parameters d a , d a , r a and r a are related as follows: d a = d a p^./p^ a ; r a = 



-r* a d a /d* a ; D a = d a d* a [H]. 

a superconductor due to the interference of electronlike and holelike 
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particles with different spin-dependent phase shifts. To demonstrate 
this, let us consider diagrams in Fig. 3, corresponding to Andreev 
reflection of an electron with the spin projection a and the energy 
less than |A| transmitted from a ferromagnet into a superconductor. 
The amplitude a(e, 6 a ) is: 

eh \-i i ~* ~ aeh ohe 



a(e, 9 a ) = d a d*_ a (3^_ a [l + r*_ a r a p^_ a P. 



-a, a 



,/n n V F / V F iP r a geh 



>eh /Q he 
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The corresponding probability of Andreev reflection is: 
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It follows from formula (TTTj) that in the absence of the interferential 
term Q the probability of Andreev reflection is a constant (indepen- 
dent of the energy e) quantity. The interference of electronlike and 
holelike particles reflected by the pair potential and interface results 
in the formation of Andreev surface bound states. At 9 a =0 the 
maximum in the probability of Andreev reflection is at e = ± |A| 
[T9] . At 6 a = ±7r/2 Andreev surface bound states with the width 
T equal to: 

r (l-v^fRi)|A| 



(19) 
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are formed at e = on the Fermi level. The peak in the differential 
conductance of an FIS contact at the zero voltage may be used 
to determine the polarization of strong ferromagnets by comparing 
experimental data with those calculated according to formula (TTU|) . 

Thus, in the present paper the ballistic conductance of the point 
FIS contact is calculated. The dependence of Andreev surface bound 
states on the spin-dependent phase shifts of the electron states re- 
flected from and transmitted through the potential barrier is found 
for the interface with finite transmission. By the example of a rect- 
angular potential barrier it is shown that these states are manifested 
in the peaks of the dependence of the conductance of the FIS contact 
on the applied voltage. 
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